Colossal anisotropy in diluted magnetic topological insulators 
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We consider dilute magnetic doping in the surface of a three dimensional topological insulator 
where a two dimensional Dirac electron gas resides. We find that exchange coupling between mag- 
netic atoms and the Dirac electrons has a strong and peculiar effect on both. First, the exchange- 
induced single ion magnetic anisotropy is very large and favors off-plane orientation. In the case of 
ferromagnetically ordered phase we find a colossal magnetic anisotropy energy, of the order of the 
critical temperature. Second, a persistent electronic current circulates around the magnetic atom 
and, in the case of a ferromagnetic phase, around the edges of the surface. 

PACS numbers: 



The prediction [lj-Ql and subsequent experimental 
confirmation^ Q of topological insulators opens an en- 
tirely new research field with fascinating possibilities 0- 
9]. Three dimensional topological insulators, like Bi2Te3 
and Bi2Se3 , host a helical two dimensional electron gas 
in their surfaces governed by the Dirac equation in which 
the spin of the electrons is strongly coupled to their mo- 
mentum, in contrast with graphene, for which the com- 
ponents of the Dirac spinor are associated to the atomic 
sublattice. In this way, spin-orbit coupling is the dom- 
inant term in the Hamiltonian and makes this system 
different from most of electron gases studied so far. 

The helical electron gas has exotic and non-trivial re- 
sponses to external perturbations, like applied electric 
and magnetic fields, the proximity of superconductors [loj 
and ferromagnets If[ 12[ . Diluted Magnetic doping is 
particularly ap pea lin g 13l420| | and has been demonstrated 
experimentally [lj, |15|. On one side, it could result in 
bulk ferromagnetic order coexisting with the topological 



order[I6t |20| |. very much like conventional diluted mag- 
netic semiconductors preserve the electronic properties 
of the parent compounds 2l|-|24(. On the other side, the 
helical electron gas is expected to affect strongly the spin 
dynamics of the magnetic atoms. In particular, because 
of the large spin-orbit coupling, the helical electron gas 
might induce a large single ion magnetic anisotropy on 
the magnetic atoms, as it happens in the case of Co in 
Pt [25j]. This diluted magnetic system could be gated [2^], 
which could permit to control both the properties of a col- 
lective ferromagnetic phase 12711 as well as the properties 
of a single magnetic atom [28|, and it could be probed 
with STM, which permits additional insight and control. 

Here we study how the spin of magnetic atoms and the 
Dirac electron gas in the surface of a topological insula- 
tor influence each other. For that matter, we assume that 
the surface of the topological insulator is weakly doped 
with magnetic atoms (see figure 1). We assume that the 
Fermi energy is not changed. This could be achieved ei- 
ther with isoelectronic magnetic dopants, by controlling 



the charge with a gate[27[ or by additional compensat- 
ing non-magnetic dopants [2^]. We find that, due to its 
coupling to the Dirac electron gas, a single magnetic ion 
has an additional source of magnetic anisotropy which 
favors off-plane alignment. We study the case of a den- 
sity of many magnetic dopants ferromagnetically ordered 
due to indirect exchange coupled mediated by the Dirac 
electron gas. We find that the magnetic anisotropy en- 
ergy density of this phase is the same than the energy 
gained by the system due to magnetic order. We refer 
to this as colossal magnetic anisotropy, in analogy to a 
similar result obtained for Pt chains [30]. Finally, we 
show that magnetic dopants induce charge currents in 
the Dirac electron gas. 

We consider the following first quantization Hamilto- 
nian for local spins f2/ and Dirac electrons: 



H = v F T-p + ■ yVflf- fr|)f2j 



(1) 



In the first term, the Dirac Hamiltonian for two dimen- 
sional electrons, the electron momentum p = (p x ,Py) lies 
in the plane. In the second term, the spin S = \t 
is exch ang ed coupled to O/, which are classical unit 
vectors [17j representing the spin orientation of the mag- 
netic atom located in 77 and J(r) is an exchange po- 
tential of range a, the interatomic distance. We assume 
half- filling, i.e., only the low energy band is occupied. 

In the following we consider two situations: the sin- 
gle spin case and the many spin ferromagnetic phase in 
the virtual crystal approximation [2l|. We treat the sin- 
gle spin case numerically. For that matter, we represent 
the Hamiltonian ([I]), for a single spin f2, in the basis of 
eigenstates of the unperturbed Dirac Hamiltonian: 



Ak-r 



'2A 



(2) 



where A is the sample area, k — \k\(cos<f>k, siruf)k), 
s = ±1. The corresponding unperturbed energies read 
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FIG. 1: . Color online. Scheme of the surface of a topological 
insulator weakly doped with magnetic atoms, without mag- 
netic order (top panel) and in the ferromagnetic pase (bottom 
panel) 



e s (&;) = shvF\k\. In this basis, the Hamiltonian of a Dirac 
electron gas interacting with a single spin reads: 

' /(/,•■■/,' !<>• ( I' 

s.k 



,(k)6 a 



s'h,P + ^(k-k')n.(Fl^ s , rkf ) (3) 



where J(k) is the Fourier transform of J{r). We take 



J(r) 



<r/ a y 



where j has dimensions of energy 



times area, so that J(k — k') 



J_ p -(\k-k'\a/2) 2 



We 



consider N ~ 900 points in a square of area A cen- 
tered around k — (0,0). Here A is a momentum cutoff 
which satisfies Aa ~ 1. For a given spin orientation O 
we obtain the 2N eigenstates E n (Q) of eq. At half 
filling the multielectron ground state energy is given by 
Ec(&)) = J2n=i ^tP) • As expected by symmetry, we 
find that Eq does not depend on the in-plane components 
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FIG. 2: Prefactor 7(aA), for the magnetic anisotropy energy 
for a single magnetic atom, as a function of Aa, calculated 
numerically, for different values of A. The points collapse 
into a single curve 



perpendicular to the surface. In that situation, the in- 
direct exchange coupling is ferromagnetic [lf|, so that 
we can expect them to order ferromagnetically below a 
certain critical temperature that we estimate below. We 
study this ferromagnetic phase in the virtual crystal ap- 
proximation (VCA) [21 - 24 1 , which is equivalent to replace 
the inhomogeneous and non-collinear exchange potential 
in eq. ([1]) by an constant exchange field, proportional to 
the concentration of magnetic atoms, c = N- lm p/ 'A, whose 
spins all point along the same direction Cl. The effective 
Hamiltonian is now translational invariant and, in the 
momentum space, reads 



A - 

%vca = hv F k ■ f + —VL ■ 



(5) 



where A = jc. The off-plane f2 z component plays 

a role very different from Q x , fl y . The latter can be 

gauged away, by choosing a constant vector potential 
n x>y , being a function of ft z only. From our numerical Aj> = with b = ^ y (see figure 3 )_ Since this 

results we infer that the single ion magnetic anisotropy 

is given by: 



8£ = E G (n z )-E G (x) = - 1 (Aa)^—A 3 n 2 z (4) 

nvF 

where 7(Aa) is shown in figure ©. This is one of the 
main results of this manuscript: because of its coupling 
to the Helical electron gas, a single magnetic atom ac- 
quires a off-plane magnetic anisotropy. We estimate the 
single ion magnetic anisotropy taking j — J a 2 , Aa ~ 1, 
and, following reference 16], J = 0.5eV and HvfA = 0.1 
eV. We obtain 5£ ~2.5 meV. This sin gle ion term com- 
petes with RKKY interactions IB, IM 13 an d with the 
intrinsic single ion anisotropy that arises from the inter- 
play of the d electron spin orbit and their crystal field 
splitting. The single ion anisotropy measured in inelastic 



2v f 

vector potential is constant, it has no effect on the phys- 
ical properties of the system. In particular, the total 
energy does not depend on £l Xl y, as it happens in the 
single ion case. In contrast, the perpendicular magneti- 
zation can not be gauged away and it has an important 
effect on the energy bands: it opens a gap of magni- 
tude A (see figure 3). The new energy bands read now 



e s (k, H z ) = sy jA 2 Q z + (hvFk) 2 . We calculate the mul- 
tielectron ground state energy, with respect to the non- 
magnetic ground state. To leading order in the cutoff A 
we obtain 
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(6) 



scanning tunneling spectroscopy [31 . 

We now consider the case of many magnetic atoms. 
Due to the single ion anisotropy, we expect them to point 



This equation can be interpreted in two different ways. 
On one side, this is the energy that the system gains in 
the ferromagnetic phase, with perpendicular easy axis, 
with respect to the non magnetic phase. The minus sign 



denotes that the helical gas is a paramagnet, in spite of 
the fact that the density of states vanishes at the Fermi 
energy. Equation ^ permits to derive an expression for 
the transition temperature. For that matter, we expand 
the free energy of the magnetic atoms around Q z = . 
To leading order, we find an entropic term, that favors 
the non-magnetic phase, and the electronic term, which 
favors magnetic order [24| : 
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where xo — 3gjr 1 ^ ^ e cr itical temperature the pref- 
actor of SI 2 changes sign: 



S(S + l)fc 

kbTc = A 

bnvF 



(8) 



On the other hand, equation (|6]) can be interpreted as 
the magnetic anisotropy energy associated to a collective 
rotation of the magnetization from in plane to off-plane. 
Thus, we can write: 



5 A- 

A S(S 
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(9) 



This result shows the peculiar nature of a ferromagnetic 
phase in the surface of a diluted magnetic topological in- 
sulator. The magnetic anisotropy and the Curie temper- 
ature are the same. In other words, a rotation of the col- 
lective magnetization from its off-plane easy axis is as ex- 
pensive energetically as a transition to the non-magnetic 
case. Thus, a collinear magnetic phase is only possible 
when the magnetic moment is perpendicular to the sur- 
face. We refer to this as colossal magnetic anisotropy, 
in analogy with the magnetic phase predicted for Pt 
nanowires '30]. Equation ([9]) points out how radically 
different is the ferromagnetic phase in a diluted magnetic 
topological insulator compared to from the ferromagnetic 
order in metals and diluted magnetic semiconductors, 
where the Curie temperature and magnetic anisotropy 
are very different [22[, |2J|. This result arises from the 
fact that the band energy and the spin orbit coupling are 
the same in the Helical electron gas. 

We now estimate the Curie temperature for a di- 
luted magnetic topological insulator. For that matter 
we consider Bi2Te3, doped with iron, for which 5 = 2. 
We assume j — Ja 2 so that we can write kbT c = 
S(S g r 2f ^ Aa2 - We assume Aa = 1 and hv F /a = 0.1 eV. 
For J = 0.5eV and ca 2 = 0.01, we obtain k b T = 25meV, 
i.e., room temperature. Notice that this prediction scales 
with J 2 , so that a 10-fold reduction of J results in a 100- 
fold reduction of T c . In addition, a mean field prediction 
of T c is likely to overestimate the real transition temper- 
ature. 

We now consider the effect of the magnetic moments on 
the Helical electron gas. It has been shown that a single 
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FIG. 3: Color on-line. Energy bands in the non-magnetic 
phase (gapless bands) and in the magnetic phase with O = z 
(gapped case) and Q = x (shifted gapless bands) 



magnetic impurity in the surface of a topological insu- 
lator opens a gap in the local density of states 16, HH 
and, in the case of the magnetic spin pointing perpen- 
dicular to the surface, it generates a circulating in-plane 
spin density around the magnetic atom. Importantly, the 
spin and the current operators are the same, modulo a 
constant, in the case of Dirac electrons: 



j — — -Hip — eA) = —evpd 
SA 



(10) 



In consequence, the in- plane circulating spin density ob- 
tained in references ( la. Il8|) can be also interpreted as a 
circulating charge current around the magnetic impurity. 
In turn, the circulating charge current bears an orbital 
magnetic moment perpendicular to the plane. Thus the 
off-plane spin induces an off-plane orbital moment. If two 
of such magnetic atoms are close to each other and co- 
polarized, the sum of their charge currents is such that 
there is no current flow in between them, but there is 
current flowing around them. The argument can be ex- 
tended to the entire surface. Thus, we expect that in 
the edges of the surface of a ferromagnetically ordered 
diluted magnetic topological insulator there is a circu- 
lating charge current in the ground state. This persis- 
tent current induced by a space dependent magnetic or- 
der parameter is similar to the one recently predicted for 
graphene with magnetic order in the edges [33j. 

We now provide a general argument to show that, in 
the presence of a off-plane varying spin density, the Dirac 
electrons must have an in-plane charge current. We con- 
sider the Dirac electrons in the presence of an exchange 
potential V = iQ, z {x)r z . Thus, the components of the 
eigenstates ^f(x,y) can be written as e lkyV 4> a (x). It can 
be easily seen that, for a state (x, y) with energy E 
and momentum k y the following identity is satisfied, for 
an arbitrary £l z (x): 



2ES y (x) + hvpk y n(x) = hvpd x S z {x) 



(11) 



where S a {x) = jEy^W^^W and n ( x ) = 
4>*(x)(f> a (x) are the spin and charge densities associ- 
ated to { x iU)i respectively. Since a given state with 
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energy E and momentum k y has a degenerate partner 
with momentum — k y , when summing over all the states 
with the same energy E, and using equation (fTDj) . we are 
left with the equation: 

{Jy{x)) E = ^{d x S z {x)) E (12) 

where () e stands for sum over all the states with energy 
E. Thus, as we move from outside the sample x < to 
the magnetic sample x > 0, there is a finite (d x S z (x))E 
which implies current flow perpendicular to the sample 
boundary. A similar analysis accounts for the current 
vortices induced by a single impurity. 

In conclusion, we study a three dimensional topologi- 
cal insulator with diluted magnetic atoms in its surface. 
We predict that magnetic dopants exchanged coupled to 
the Dirac electron gas in the surface of a topological insu- 
lator acquire a magnetic anisotropy that favors off-plane 
orientation. We find that the magnetic anisotropy en- 
ergy and the energy gain due to magnetic order are the 
same. As a result of this colossal magnetic anisotropy, 
the ferromagnetic phase can only exist when the spins 
are oriented perpendicular to the plane. We find that 
magnetic dopants induce in-plane charge currents. In the 
ferromagnetic phase, charge current flows in the edges of 
the sample. 
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